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I. INTRODUCTION 

The transition from the confinement to the decon- 
finement phase is one of the most interesting fea- 
tures of QCD at finite temperature. This transition 
separates a low-temperature confinement phase from 
a high-temperature - quark-gluon plasma - phase, 
where color charges should be deconfined. The ex- 
istence of this transition has been confirmed by recent 
observations of the collective effects in ultrarelativistic 
heavy- ion collisions (see, e.g., the review 

The non-perturbative study of gauge variant prop- 
agators is of interest for various reasons. These prop- 
agators are expected to show different behavior in 
each phase and, therefore, to serve as a useful 'or- 
der parameters', detecting the phase transition point 
T c . One expects that their study can shed the light 
on the mechanism of the confinement-deconfinement 
transition [2-5]. Another reason is that for the reli- 
able phcnomcnological analysis of high-energy heavy- 
ion collision data, it is important to obtain informa- 
tion on the momentum dependence of the gluon prop- 
agators, especially in the infrared region (see, e.g., p- 
fLOl]). Finally, the non-perturbatively calculated lattice 
propagators are to be used to check the correctness 
of various analytical methods in QCD, for example, 
Dyson-Schwinger equations method at finite temper- 
ature dMl. 

Last years a number of papers have been dedicated 
to the lattice study of the finite temperature SU(2) 
pure gauge gluon p ropagators in Landau gauge (see, 
e.g., papers [13M20I ]). However, the finite temperature 
SU(3) gluon propagators - especially in the presence 
of dynamical fermions - are much less studied (see 
papers [H, [21H23) for the pure gauge theory and [24| 
for the full QCD). 

In this paper we study the transverse (magnetic) 



and longitudinal (electric) gluon propagators both in 
the pure gauge SU(3) theory and in the full QCD 
with the number of flavors np = 2. The main goal 
of our work is to compare the momentum behavior of 
these propagators very close to the transition point, 
slightly above T c (deconfinement phase) and slightly 
below T c (confinement phase). We compare our re- 
sults for the gluon propagators computed in gluody- 
namics and in full QCD at coinsiding values of the 
ratio T/T c where T c is a corresponding transition 
temperature. For both theories we have chosen tem- 
peratures T very close to the corresponding transition 
point : T/T c = 0.97 and T/T c = 1.02 . Of special 
interest in this study are the infrared mass scale pa- 
rameters ('screening masses'): tue ('electric') and tum 
('magnetic'). 

We apply effective gauge fixing algorithm - simu- 
lated annealing - and generate few gauge copies for 
every configuration to reduce Gribov copy effects. For 
gluodynamics we apply additionally flips between sec- 
tors of the Polyakov loops [25[ . The flips were applied 
to this theory for the first time. Also we check finite 
volume effects comparing results for two lattice sizes. 



Section ILTl contains main definitions as well as some 
details of simulations and gauge fixing procedure we 
use. Volume and temperature dependence of the prop- 
agators for both pure gauge theory and for QCD are 
discussed in Section IIHI Section IIVI is dedicated to 
the discussion of the screening masses and Section [V] 
is reserved for conclusions and discussion. 
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II. GLUON PROPAGATORS: THE 
DEFINITIONS AND SIMULATION DETAILS 

For the study of the gluon propagators in the pure 
gauge theory we employ the standard Wilson action 
Sw with j3 — 6/<7q where go is a bare coupling con- 
stant. To define the spacing a as a function of /3 
in the pure gauge case we have used Necco-Sommer 
parametrization [26| with the popular choice of the 
Sommer scale r$ — 0.5fm. In what follows we will 
refer to this theory as np = theory. 

To study the effect of the quarks on the gluon propa- 
gator we computed propagators on configurations gen- 
erated with the gauge action Sw and rip — 2 dynami- 
cal flavors of nonperturbatively 0{a) improved Wil- 
son fermions (clover fermions). The configurations 
were produced by the DIK collaboration [27J with 
BQCD code [28j. The dimensionless quantities ro/a 
and m T ro are taken from results of QCDSF collabora- 
tion. To convert them into physical units we use the 
value r = 0.467(15) fm obtained by QCDSF gg|. The 
improvement coefficient csw was determined nonper- 
turbatively [3(| • In the rest of the paper this theory 
will be refered to as np = 2 theory. 

Our calculations were performed on the asymmet- 
ric lattices with lattice volume V — L4, ■ L\ , where 
Li is the number of sites in the 4th direction. The 
temperature T is given by T = 1 / aL^ where a is the 
lattice spacing. We employ the standard definition of 
the lattice gauge vector potential A xfjL [3l| : 
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where link variables U x>1 G SU (3) transform under 
gauge transformations g x £ SU(3) as follows : U xfl & 

The Landau gauge condition is given by (dA) x = 

2^=1 (A>m - Ax-fop) = which is equivalent to 
finding an extremum of the gauge functional 



M9)= I ^> 

XfJ, 



(2) 



with respect to gauge transformations g x . 

The bare gluon propagator D?^(p) is given by 



Df v {jp) = ^Ul{k)A b v {-k) 
9o x 



(3) 



where A{k) represents the Fourier transform of the 
gauge potentials according to Eq. (TTJ) after having 
fixed the gauge. The physical momenta p are given by 
Pi = (2/a) sin (jrki/L), p^ — (2/a) sin (irk^/L^), fcj G 
(-L/2,L/2],k 4 e(-L A /2,L 4 /2}. 



In what follows we consider only soft modes p^ = 0. 
The hard modes (p 4 7^ 0) have an effective thermal 
mass 2nTn and behave like massive particles. 

On the asymmetric lattice there are two tensor 
structures for the gluon propagator 



32] 



DfM = 5 ab (PT(P)D T (P) + P^(p)Dl( P )) , (4) 

where (symmetric) orthogonal projectors P^i (p) are 
defined at p = (p ^ 0; P4 — 0) as follows 



pL(p) = 1 ; i$(p) = o 



(5) 
(0) 



Therefore, two scalar propagators - longitudinal 
D L (p) and transverse D T (p) - are given by 



8 3 



a—l i—l a—l 

(7) 

For p — propagators D T (0) and D L (0) are given by 

a—l i—l a—l 

(8) 

The transverse propagator D T (p) is associated to 
magnetic sector, and the longitudinal one D L (p) - to 
electric sector. 

In the case of the np = theory we employ for 
gauge fixing the Z(3) flip operation as has been pro- 
posed in [25j . It consists in flipping all link variables 
U xfl attached and orthogonal to a 3d plane by multi- 
plying them with exp {±27ri/3}. Such global flips are 
equivalent to non-periodic gauge transformations and 
represent an exact symmetry of the pure gauge action 
Sw- At finite temperature we apply flips only to di- 
rections fi = 1, 2, 3. As for the 4th direction, we stick 
to the sector with \arg P\ < 7r/3 which provides max- 
imal values of the functional F. Therefore, the flip 
operations combine for each lattice field configuration 
the 3 3 distinct gauge orbits of strictly periodic gauge 
transformations into one larger gauge orbit. 

All details of our gauge fixing procedure - FSA algo- 
rithm - are described in our recent papers [l9l . |33T - [35| . 
The features of the simulated annealing algorithm ap- 
plication specific for SU(3) group are essentially same 
as m 36] where this algorithm was applied to SU (3) 
theory for the first time. For every configuration we 
produce two gauge copies per flip-sector; therefore we 



have in total N c , 



54 copies. We take the copy 



with maximal value of the functional ([2]) as our best 
estimator of the global maximum and denote it as best 
("6c ") copy. 
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In the case of the np — 2 theory the global Z(3) 
transformations do not anymore make part of the 
symmetry group of the action (only one flip-sector re- 
mains). In this case we have made N copy — 10 gauge 
copies. 

To suppress 'geometrical' lattice artifacts, we have 
applied the "a-cut" [37}, i.e. irki/L s < a , for every 
component, in order to keep close to a linear behav- 
ior of the lattice momenta pi ~ (2nki)/(aL s ), ki G 
(—L s /2,L s /2\. We have chosen a = 0.5. Obviously, 
this cut influences large momenta only. We did not 
employ the cylinder cut in this work. 

The values of T/T c , (3, k and lattice sizes are given 
in Table [I] Number of independent configurations in 
all cases equals ~ 200. 





n F = 


np = 2 


T/T c 


1/4 • Z/ 3 


P 






0.97 


8- 16 3 


6.044 


5.25 


0.1341 


1.02 


8- 16 3 


6.075 


5.25 


0.1345 


0.97 


8 • 24 3 


6.044 


5.25 


0.1341 


1.02 


8 • 24 3 


6.075 


5.25 


0.1345 



TABLE I: Run parameters. 



For np = theory the critical value /3 C (correspond- 
ing to the 1st order phase transition) has been taken 
from gl]. 

In the case of np — 2 theory the simulations were 
made [13] at fixed (3 — 5.25 and varying k. The 
pseudocritical values k c (and thus the transition tem- 
peratures T c ) of the deconfining and chiral symmetry 
restoration transitions have been determined [27| by 
the maximum of the Polyakov loop and chiral conden- 
sate susceptibilities, respectively. It has been found 
[27| that these two maxima coinside within numerical 
precision at k c — 0.1343. 

The pion is heavy at values of k we have chosen 
(m-n ~ 1 Gev). However, it is known that the effect 
of the quarks is rather strong even at these values of 
TYL-K . The transition is a crossover, rather than a 1st 
order phase transition. Furthermore, the transition 
temperature is shifted to substantially lower values in 
comparison to the gluodynamics [27J, |39|, |40j and the 
string breaking phenomenon is observed at T < T c 

Mm- 



III. GLUON PROPAGATORS : NUMERICAL 
RESULTS 

We define the renormalized propagators D^ en (p) 
and D^ en (p) in such a way that their dressing func- 
tions are equal to unity at the normalization point 
fi = 2.5 Gev. In the rest of this paper we will omit 



the subscript Yen'; therefore, Dt{p) and D^(p) will 
denote transverse and longitudinal renormalized prop- 
agators, respectively. 





n F — 


np — 2 


T/T c 


Z/4 ■ Z/ 3 


Zm 


z E 


Zm 


Z E 


0.97 
1.02 


8- 16 3 
8- 16 3 


1.03(1) 
1.01(1) 


1.00(2) 
1.00(1) 


1.12(2) 
1.13(2) 


1.14(2) 
1.16(2) 


0.97 
1.02 


8 • 24 3 
8 ■ 24 3 


1.04(1) 
1.04(1) 


1.00(1) 
1.03(1) 


1.14(1) 
1.15(1) 


1.14(1) 
1.17(1) 



TABLE II: Renormalization constants Z M /p. 



To compute the renormalization constants Zm and 
Zp we fitted the propagators at momenta p in the 
range between 2.3 and 2.7 GeV to a polinomial func- 
tion. They are shown in Table [TTI for both np = 
and np = 1 theories. As expected, there is no volume 
dependence of these constants. In the case of the pure 
gauge theory all Zm and Zp are very close to unity, 
which is in agreement with the results for the bare 
gluon propagators at zero temperature in pure gauge 
SU(3) theory 43]. In the np = 2 case all renormaliza- 
tion constants are substantially larger. This deviation 
(~ 15%) might be explained, at least, partially, by 
the quark contribution. Another possible explanation 
might be the larger spacing (or smaller j3) effect. 

For both theories every constant, i.e., Zm and Zp, 
is practically the same above and below T c . This 
means that ultraviolet parts of Dt(p) and D L (p) are 
not 'phase-sensitive'. 

A. On Gribov copy effects 

In our study of the Gribov copy effects we follow the 
approach described briefly in Section [Til an d - in more 
details - in our papers [l! [H HlMlill dedicated 
to the Gribov copy effects in SU(2) gluodynamics at 
zero and nozero temperatures. 

We find that in the case of the nj? = theory, simi- 
lar to SU (2) gluodynamics, the Gribov copy influence 
still remains a serious problem in the lattice calcula- 
tions, at least, at small momenta. These effects arc 
very strong for Dt{p) in the infrared and weak for 
Z?l(p) for all momenta. 

It is worthwhile to note that the effect of different 
flip-sectors plays a very important role in this case, 
even more important than additional gauge copies in 
each sector. The details of this study will be published 
elsewhere. 

In contrast, in the case of dynamical fermions (np = 
2) the Gribov copy effects look completely negligible 
as compared to the pure gauge case. The propaga- 
tors do not change within errorbars ("Gribov noise"). 
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However, it is worthwhile to note that our lattice sizes 
are comparatively small. One cannot exclude that 
these effects might become stronger for larger lattice 
sizes and/or smaller lattice spacings. 



B. Volume dependence 

We find that for the transverse propagators Dt(p) 
the finite volume effects are rather strong at zero mo- 
mentum but they become very weak with increasing 
\p\. As an illustration, in Figure Q] we show the mo- 
mentum dependence of the np — and np = 2 trans- 
verse propagators at comparatively small momenta 
(\p\ < 1.5 GeV) at T/T c = 1.02 on the lattices 8 ■ 16 3 
and 8 • 24 3 . No volume dependence is observed also 
for \p\ > 1.5 Gev. The similar situation takes place 
for T/T c = 0.97. 

The volume dependence of the longitudinal propa- 
gators Dl(p) is different. One can see from Figure [2] 
that the sign of the finite volume effects differs from 
that for Dt{p) case : while the transverse propagator 
decreases with increasing volume, the longitudinal one 
is increasing. This is in agreement with SU (2) theory 
[l9j |. Another observation is that the finite volume 
effects are much more substantial for np = theory 
than for np = 2 one. At p — the quenched the- 
ory propagator increases by factor two when L s in- 
creases from 16 to 24 as can be seen from Figure [5] at 
T/T c = 0.97. This increase is similar to that observed 
in SU(2) theory [II 01]. The finite volume effects de- 
crease fast with increasing \p\ both for np — and for 
np = 2 theories. At T/T c = 1.02 the corresponding 
figure looks in a similar way. 

In what follows we show data for L s = 24 lattices 
only. 
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FIG. 1: The momentum dependence of the transverse 
propagator D T {p) at T/T c = 1.02 for \p\ < 1.5 Gev for 
two lattice sizes: L s = 16 and 24. 
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C. Temperature dependence 

We find that the transverse propagators Dt(p) are 
not very much sensitive to the crossing of the transi- 
tion point T c , at least, at p ^ 0. In Figure [3] where 
we present results for Dt{p) at T/T c = 0.97 and 
T/T c = 1.02 for both n/ = and u/ = 2 theories. We 
observe that in the case of n/ = theory the temper- 
ature dependence is very weak at all momenta. (Sim- 
ilar conclusion was made in (l8|.) The same is true 
for nf — 2 theory at nonzero momenta. But, contrary 
to rif = case, for p = substantial change in the 
propagator is observed as can be seen from Figure [3J 
This is one of the clear effects of dynamical quarks. 

Another interesting observation is that one can 
not see any substantial influence of the dynamical 
fcrmions on the momentum dependence of the trans- 
verse propagators for \p\ ^ 0. Moreover, for the 
momenta \p\ ^ 0.8 Gev all four propagators coinside 



FIG. 2: The momentum dependence of the longitudinal 
propagator D L (p) at T/T c = 0.97 for |p| < 1.5 Gev. 



within errorbars. It would be interesting to see how 
this situation will change for smaller values of the pion 
mass m T . The values of propagators obtained at p = 
suggest that in nf — 2 theory the transverse propaga- 
tor is infrared enhanced in comparison with the one 
in nf =0 theory. This is in contrast with T = re- 
sults [H, Ell ■ But this conclusion should be checked 
in future studies on larger volumes. 

In contrast, the longitudinal propagators 
demonstrate much more pronounced temperature de- 
pendence for both theories, especially in the infrared 
region |p| < 1.5 Gev. For both theories, i.e., np = 
and np — 2, the propagators Dl(p) differ significantly 
also for nonzero values of momenta when temperature 
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FIG. 3: The momentum dependence of the transverse 
propagator Dt(p) at two temperatures for \p\ < 1.5 Gev. 



FIG. 5: The momentum dependence of the longitudinal 
propagator Dl(j>) at two temperatures for \p\ > 1.5 Gev. 
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FIG. 4: The momentum dependence of the longitudinal 
propagator Dl(p) at two temperatures for \p\ < 1.5 Gev. 



FIG. 6: The momentum dependence of Dt(p) and Dl(p) 
at two temperatures for |p| < 1.5 Gev. 



T crosses the transition point T c (see Figure |4j . Fast 
change of Dl(j>) across transition in np =0 theory 
was found recently in [is) . 

Fermions influence the longitudinal propagators up 
to values of momenta \p\ ~ 1.8 Gev what can be seen 
in Figure [S] For larger values of \p\ we see no differ- 
ence between Dl(p) computed in np = and np = 2 
theories as well as no temperature dependence. 

Finally, we compare the momentum dependence of 
the transverse and longitudinal propagators for both 
theories at different temperatures. In Figure [S] the 
comparison of Dt (p) and D l (p) in the infrared region 
(\p\ < 1.5 Gev) for np — 2 theory is present. 

Evidently, there is a large difference between Dt(p) 
and Dl(p) at small momenta. However, this differ- 
ence quickly decreases with increasing \p\ and becomes 



zero (within errorbars) for |p|^1.2 Gev. 



IV. ON THE SCREENING MASSES 

One of the interesting features of the finite temper- 
ature physics is the appearance of the infrared mass 
scale parameters : mp ('electric') and uim ('mag- 
netic'). These parameters (or 'screening masses') de- 
fine screening of electric and magnetic fields at large 
distances and, therefore, control the infrared behavior 
of Dp(p) and Dp{p)- The electric screening mass niE 
has been computed in the leading order of perturba- 
tion theory long ago: m 2 E /T 2 = ( 1 + n F /6)g 2 {T) 
for 517(3) theory with n F flavors [UlU]. But at the 
next order the problem of the infrared divergencies has 
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FIG. 7: The momentum dependence of the inverse propa- 
gators D~[ 1 (p) for two different temperatures. Solid lines 
correspond to fit to eq.{TU]). 



been found. On the other hand, the magnetic mass 
rriM is entirely nonperturbative in nature. Thus a 
first-principles nonperturbative calculations in lattice 
QCD should play an important role in the determina- 
tion of these quantities. 

In the deep infrared region the momentum depen- 
dence of the longitudinal propagator is expected to fit 
the pole-type behavior 



D L (p) 



c 



- jr 



p ~ 



(9) 



Indeed, in the case of the pure gauge SU (2) theory we 
observed the linear behavior of the inverse longitudi- 
nal propagator Dj^(p) as a function of p 2 for small 
enough values of momenta [19j . Moreover, this linear 
dependence has been observed both above and below 
T c and even at T ~ T c . 

In the present work our lattice sizes L s do not per- 
mit us to reach small enough values of p 2 where the 
linear p 2 -behavior of (p) could be found. In- 
stead we have used for the fit a somewhat symplified 
Stingl-like parametrization of the longitudinal propa- 
gator m i5(j 



DZ 1 (p) = C- 1 -(m 2 E +p 2 + b\p 



(10) 



We emplyed this formula in the comparatively large 
momentum interval, up to p 2 ^2 Gev 2 , where 6 data 
points were used for the fit for every data set (see Fig- 
ure [7]) . The parameter b indicates the deviation from 
the simple pole- type behavior. We use the parameter 
ffiE as an estimator of the electric screening mass mg. 

The values of the fit parameters are presented in 
the Table IIII1 In the case of the pure gauge theory 



T/T c 


riF 


rhs/T 


C 


bT 2 


0.97 





1.725(126) 


5.03(49) 


0.091(14) 


0.97 


2 


2.09(15) 


3.17(14) 


0.054(6) 


1.02 





1.851(77) 


3.24(15) 


0.050(5) 


1.02 


2 


2.24(9) 


2.68(16) 


0.045(7) 



TABLE III: Values of parameters the, C and b obtained 
from fits to eq.fTOj). 



we can compare our value of m^/T calculated at 
T/T c = 1.02 with m E /T obtained in [H[ for nearby 
temperature T/T c — 1.05 . Both values are well con- 
sistent within errorbars. Instead, in the unquenched 
QCD case our value of m E /T at T/T c = 1.02 differs 
by factor 1.5 approximately from mp/T calculated 
in np =2 QCD with KS fermions at the same value 
of T/T c [11]. The origin of this deviation is still to to 
be clarified. 

As is well-known, one can hardly rely on the pertur- 
bation theory at T ~ T c . However, it is interesting to 
note that the ratio of two values m 2 E /T 2 (for np = 2 
and rip = 0) is consistent with the factor (1 + np/Q). 

Our main observation is that the dimensionless pa- 
rameter rhp/T is practically not sensitive (within 
errorbars) to the crossing of the transition point T c . 
This statemant is true not only for np = 2 theory 
where the crossover (or higher order phase transi- 
tion^ expected but also for the pure gauge theory 
where the existence of the 1st order phase transition 
is firmly established. Therefore, we conclude that the 
electric mass can hardly be considered as an 'order 
parameter' indicating the transition point T c . 

Note, that a similar situation have been observed in 
the case of the finite temperature pure gauge SU{2) 
theory where the values of the parameter m^/T 
coinside for T/T c = 0.9 and T/T c ~ 1 [lg. 

In contrast, dimensionless factors C are much more 
sensitive to the crossing of the transition point T c . 
This is true for both np — and np ~ 2 theories. 
The difference between np = and np = 2 cases is 
more quantitative than qualitative : the temperature 
variation of the parameter C in the pure gauge theory 
is much more strong as compared to the full QCD 
case. 

In fact, the only qunatity which demonstrates clear 
difference between np = and np — 2 cases is the 
dimensionless parameter bT 2 . Indeed, for the pure 
gauge theory the temperature jump of bT 2 is rather 
strong, while for the np = 2 case values of bT 2 above 
and below T c coinside within errorbars. 

The calculation of the magnetic screening mass nriM 
looks a somewhat more delicate problem. In the case 
of pure gauge SU(2) theory it has been shown that at 
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T > T c the transverse pr opag ator Dt{p) has a max- 
imum at p ^ [TH [H, Il8ll This has been found 
also in SU(3) gluodynamics [18| . Moreover, in our re- 
cent paper [19] it has been shown that Dt(p) has a 
maximum not only at high temperatures but even at 
T < T c . The position of this maximum at T close to 
T c was found to be at about 400 Mev. Therefore, the 
transverse propagator Dt(p) has a form which is not 
compatible with the simple pole-type behavior, so for 
tum another, different from pole mass, definition is 
necessary. This definition has been proposed in [l9| . 

In paper [5l| it has been suggested that the proxim- 
ity of the Gribov horizon at finite temperature forces 
the transverse gluon propagator Dt(p,Pa = 0) to 
vanish at zero three-momentum. If this is the case, 
then the finite-temperature analog of Gribov formula 
\p\ 2 /(\^ 4 + Mf { ) = l/(\p\ 2 + ml ff (f)) suggests that 
the effective magnetic screening mass m e ff(p) be- 
comes infinite in the infrared (so called, magnetic glu- 
ons 'confinement'). 

In the case of SU(3) gauge group, both for n/ =0 
and rif = 2 theories considered in this paper, we do 
not see a maximum of Dt(j>) at p ^ 0. The transverse 
propagators Dt(p) look somewhat similar to the lon- 
gitudinal ones. The most probable reason for this is 
that our volumes are not large enough. Note that 
also in the case of SU(2) theory we did not observe 
the maximum at p ^ on smaller lattices; this maxi- 
mum appeared only with increasing of lattice size L s 
(l9j . To fit Dt (p) we employed the same formula given 
in eq. (ITUl) and same fitting range as for Dl(p). We 
found that the fit is good when p = is excluded. 
Otherwise the x 2 /ndf increases up to 6 and becomes 
significantly larger than in the case of Dl(p). We con- 
sider this as an indication of the maximum to be seen 
when the volume is increased. 

We conclude that for reliable definition of the mag- 
netic screening mass uim one needs larger values of 
the lattice size L,. 



V. CONCLUSIONS 

In this work we studied numerically the behavior of 
the Landau gauge transverse and longitudinal gluon 
propagators in the pure gauge SU(Z) theory and in 
the theory with dynamical fermions with the number 
of flavors np — 2. We compare these two theories in 
the close vicinity of the transition temperature T c , 
slightly above and slightly below corresponding tran- 
sition temperature. For both theories we have chosen 
temperatures T in such a way that ratios T/T c are 
equal to 0.97 and 1.02. To our knowledge, this is the 
first study where two theories are compared in con- 
finement and deconfinement phases. 

Let us summarize our findings. 

In the case of the pure gauge theory the Gribov 



copy effects are rather strong for Dt (j>) in the infrared 
and weak for Dp(p) for all momenta. In contrast, in 
the case of dynamical fermions the Gribov copy ef- 
fects look completely negligible (the so called " Gribov 
noise"). 

The renormalization constants Zm and Zp are 
practically the same above and below T c for both the- 
ories. It confirms that ultraviolet parts of Dt{p) and 
Dl(jj) are not 'phase-sensitive'. 

Both below and above T c the finite volume effects 
are rather strong for the transverse propagators Dt(p) 
at zero momentum but they become weak fast with 
increasing \p\. This weakening is slower for n/ = 2 
theory. The volume dependence of the longitudinal 
propagators Dl(p) is different. First, the sign of the 
effect is different from that for Dt{p) case. Second, it 
is much more substantial for np — theory than for 
np = 2 one. At p = the quenched theory propagator 
increases by factor two with increasing volume as can 
be seen from Figure [5] This increase is similar to that 
observed in SU{2) theory [13, [Hf ■ However, at p ^ 
this dependence is not very much pronounced, both 
for np — and for np — 2 theories. 

We computed the electric mass mp defined in 
cq. fTU)) . It is practically not sensitive to the crossing 
of the transition point T c . This is true not only for 
np — 2 theory where the crossover (or higher order 
phase transition) is expected but also for the np = 
case where the existence of the 1st order phase tran- 
sition is well established. Therefore, mp can hardly 
be considered as an 'order parameter' indicating the 
transition point T c . 

In contrast, factors C defined in eq. (fit)]) are much 
more sensitive to the confinement-deconfinement tran- 
sition for both np — and np — 2 theories. However, 
the difference between np = and np — 2 cases is 
more quantitative than qualitative : the temperature 
variation of C in the pure gauge theory is much more 
strong as compared to QCD case. Note, that temper- 
ature dependence for mp and C computed in SU(2) 
gluodynamics [19[ was similar to those found here in 
np = theory. 

The only (dimensionless) parameter which demon- 
strates clear difference between np = and np = 2 
theories near T c is bT 2 . For the pure gauge theory the 
confinement-deconfinement variation of bT 2 is rather 
strong. Instead, for the full QCD case values of bT 2 
above and below T c coinside within errorbars. 

Calculation of the magnetic screening mass wlm 
is somewhat more delicate problem. Our caution is 
based on the experience of Tom study in the SU(2) 
theory. Indeed, in the case of the pure gauge SU{2) 
theory the transverse propagator Dt{p) has a local 
maximum at p ^ in the deconfinement and even in 
the confinement phases when lattice size L s is large 
enough and Gribov copies are properly handled. In 
our study of SU(3) theories we do not see a maxi- 
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mum of Dt(p) at p ^ and the transverse propaga- 
tors Dt{p) look similar to the longitudinal ones. It 
is very probable that this is because our lattice sizes 
are not large enough. We conclude that for reliable 
definition of wlm one needs larger values of L s and 
/3's. 
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